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Abstract 



Charged Higgs boson production via the gluon-bottom quark mode, gb — ► iff 111 , 
followed by its decay into a chargino and a neutralino has been investigated. The 
calculations are based on masses and couplings given by the Minimal Supersym- 
metric Standard Model (MSSM) for a specific choice of MSSM parameters. The 
signature of the signal is characterized by three hard leptons, a substantial missing 
transverse energy due to the decay of the neutralino and the chargino and three 
hard jets from the hadronic decay of the top quark. The possibility of detecting the 
signal over the Standard Model (SM) and non-SM backgrounds was studied for a set 
of tan/3 and rriA- The existence of 5-a confidence level regions for H discovery at 
integrated luminosities of 100 fb _1 and 300 fb _1 is demonstrated, which cover also 
the intermediate region 4 < tan (3 < 10 where decays to SM particles cannot be 
used for discovery. 



1 Corresponding author: Christian. Hansen@cern.ch 



1 Introduction 



The search for Higgs bosons is at the front-line of present research efforts in particle 
physics. While there is a single Higgs boson in the Standard Model (SM) [l] (the only 
SM-particle not yet discovered), the Minimal Supersymmetric (SUSY) extension of the 
Standard Model (MSSM) has five of them Oj- Three neutrals, the CP-even h° and 
H° (where m h < m H ), and the CP-odd A , and two that are charged conjugates of each 
other, iP* 1 . The detection of the charged Higgs bosons would unambiguously imply the 
existence of physics beyond the SM, since charged scalar states (like H ± ) do not belong in 
the SM. The combined LEP collaborations have set lower limits in a model independent 
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Figure 1: The ATLAS 5-cr discovery contour of the charged Higgs jlj. Below ~160 GeV 
the processes t — *> bH ± , H ± — > tv provides coverage for most tan/3. Above ~175 GeV 
the gb — *> bH ± , H ± — > tv covers the high tan/? region (tan/5 > 10) while the H ± — > tb 
channel covers the tan f3 < 4 region. In the intermediate tan f3 region the charged Higgs 
decays to SM particles are undetectable at the LHC. 



way on the mass of if ^bosons, M#± > 78.6 GeV for any value of tan (3 [5 . At tree level in 
the MSSM, all Higgs particle masses and couplings are determined by two parameters [§]. 
The conventional choice is to use the mass of the CP-odd neutral Higgs, tua, and the ratio 
of the vacuum expectation values of the Higgs doublets, tan (3. For the choice of MSSM 
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parameters considered here, the mass of the if ± does not differ considerably from that of 
the A . The discovery potential of H at the LHC has been investigated by both ATLAS 
|lj and CMS jZj collaborations. It has been established that for itia below ~160 GeV and 
for most values of tan (3 the charged Higgs can be discovered with 95% confidence level 
(C.L.) through the process t — > bH ± , H ± — > rv |H] (see figure EJ. Above the top quark 
mass (i.e. m#± > 175 GeV) the charged Higgs is produced via the gluon-bottom quark 
mode, gb — ► tH^. In this mass region the charged Higgs can be discovered with a 5-cr 
C.L. through H ± — > tb for tan (3 < 4, such low tan (3 values are however already excluded 
by LEP with the m t value used for the current investigation 2 , and for tan (3 > 15 [TQj and 
through H ± — > rz/ for tan/5 > 10 [H]. In the intermediate tan/? region (4 < tan/? < 10) 
charged Higgs decays into SM particles have been found to be undetectable at the LHC. 
This zone of undectectability is for the if ± — > tb channel partly due to the H ± tb Yukawa 
coupling that is ~ (m b tan (3 + m t cot j3) and has a minimum at tan (3 = a/ m t /mb ~ 7. 
Also the H ± — > tv channel becomes invisible at this region due to decreasing expected 
signal rate for decreasing tan/? values [O]. For heavy charged Higgs (rriH+ > m t + rrib) 
and for large tan (3 (tan (3 > 30 for tjih+ ~ 250 GeV and tan f3 > 50 for m H + ~ 500 GeV) 
it has been shown that the charged Higgs can be discovered at the LHC through the 
process gg — > tbH ± , H ± — » tb [T3]. 




Figure 2: The Feynman diagram (a) shows the production of the iH + final state and the 
diagram (b) shows an example of how H ± can decay to 3 leptons and several undetectable 
particles. The ® means either xt — > Xi + W + where W + — » + /i + or — >• jl + + ^ 
where /i + — > + /U + or — >• i> M + /i + where z>^ — > %5 + *V 

Among the new massive sparticles predicted in the MSSM are the charginos, xf 2 and 
the neutralinos, Xi 2 3 4> which are mass eigenstate mixtures of the electroweak gauginos 
and the Higgsinos. As normally done, the xi wm De assumed to be the stable Lightest 
SUSY Particle (LSP), which is correct unless there is a lighter gravitino. 

It is concluded in ^3] (based on J3]) that searching for decays into a chargino 
and a neutralino can be a viable method for discovery in parts of the intermediate 

2 This exclusion region can however change with the new world-averaged m t = 178.0 ± 4.3 GeV 
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region 4 < tan/3 < 10. The authors show that the processes gb — > tH + + c.c. with 
H ± — > X12X1 234 can be distinguished from the SM and non-SM backgrounds using as 
signature three hard leptons and substantial missing transverse energy from the neutralino 
and chargino decays and three hard jets from the hadronic top decay. In figure El an 
example is shown of how the charged Higgs decay in the chargino neutralino channel 
—> X12X1 2 3 4 ^ 3 ^ + jV can produce three leptons and a number, N, of undetectable 
particles in the final state. This paper reports on an investigation of the charged Higgs 
discovery potential for this channel in the intermediate tan (3 region using the ATLAS 
detector. The dominating SM background channels are gg — > ti and gg — > tiZ. A 
third SM background channel, gg — > ttj, was demonstrated in [H] to give a negligible 
contribution. The dominating non-SM backgrounds are gg — > tih, gg — > XX an d gg — > 
q,g. For this investigation integrated luminosities of 100 fb _1 and 300 fb -1 are assumed. 
HERWIG [TE], [T7j is used as Monte Carlo physics simulator and ATLFAST jTH] is used 
to perform a fast simulation of the ATLAS detector 3 . 

In section |2] the used MSSM parameters and the signal branching ratio (BR) are 
presented. In section 01 the cross sections for the charged Higgs production and the 
background channels are discussed. Section HI shows the expected number of events and 
in section |3] we describe the analysis. Finally the results are summarized and discussed 
in section IH1 



2 MSSM Parameter Point and Signal Branching Ra- 
tio 

For the analysis shown in this paper a point in the MSSM parameter space is chosen so 
that the branching ratio for Higgs decays into a chargino and neutralino are maximized. 
The same point was used for the CMS-analysis in Many independent input MSSM 
parameters are involved in the calculation of the —>■ chargino-neutralino rate. At 
tree-level, two parameters, uia and tan /5, completely specify the Higgs masses and cou- 
plings to SM particles For the charginos and neutralinos (the "inos") the tree- level 
masses and couplings to the charged Higgs bosons are determined by the parameters ttla, 
tan p, M2 and /1. Mi is assumed to be determined from M2 via gaugino unification; 
Mi = I tan 2 QwMi- The branching ratio of ino decays to leptons depend on the trilinear 
coupling, A T , which is here put to zero, and on the left- and right-handed soft slepton 
masses. It is assumed that for all three generations = rrij> L and that the two first 
generations are degenerate in mass, i.e. m gi R = m^ L R = mj> . It was chosen to max- 
imize the ino decays to leptons by setting mr to the lowest value allowed by the LEP 
results which is mi ~ 110 GeV [Tl]. The third generation slepton masses, rrif L R , are 
assumed to be 100 GeV above the selectron and smuon masses. This enhances the ino 
decays to leptons since the branching ratios of ino — > stau are then less than to the first 

3 Both packages are used within the ATHENA framework in the ATLAS Software Release 7.0.3. 
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Figure 3: The contours of constant branching ratio of H ± — > Xi 2X1 2 3 4 — * 3^ + N, where 
N represents undetectable final state particles and I = or /r* 1 , are shown in (a) for 
tan (3 = 10 and m A = 500 GeV in the M 2 vs. // plane and in (b) for M 2 = 210 GeV 
and n = 135 GeV in the tan/5 vs. m A plane. Both in (a) and (b) mi = 110 GeV, 
rrir R = 210 GeV, m g = 800 GeV, m g = 1 TeV and m t = 175 GeV and the trilinear coupling 
terms are set to zero. 



slepton generations. A heavy gluino mass, m g = 800 GeV, is assigned and all squark 
masses are set to 1 TeV. These values are used as input to ISAJET7.64 [20^ to calculate 
the total branching ratio of H ± — > X^ 2 Xi,2,3,4 "~ * 3( + N in the M 2 -/i parameter space 
setting tan /3 = 10 and m A = 500 GeV. According to the results shown in figure |Hk. 4 the 
branching ratio of charged Higgs decay to three leptons is higher than 0.4 by choosing 
M 2 = 210 GeV and \x = 135 GeV. The MSSM parameter point chosen is thus: 

• M 2 = 210 GeV, n = 135 GeV, m~ lR = 110 GeV, m fR = 210 GeV, m g = 800 GeV, 
rriq = 1 TeV. 

With these values as input ISAJET7.64 was used to calculate the total branching ratio of 
~ *■ Xi = 2Xi 2 34 ^ 3£ + in the tan (3-m A space. The results are shown in figure Eb- 



3 Cross Sections 

The next-to-leading order (NLO) cross section for the signal, gb — > tH ± , has been calcu- 
lated in for the tan/5 and m A values indicated with dots in figure 5 - The contours 

4 ISAWIG 1.200 (an ISAJET-HERWIG interface) was used for these calculations. 
5 In JHj the cross sections are given for (tan/3, mu) values. The change in the cross section when 
putting huh = rn,A is negligible. 
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Figure 4: The contours in (a) shows the NLO cross section for the signal (gb — > tH + , iH~). 
The cross sections for the (tan f3, w,a) values indicated with dots are taken from ^H] and 
in between the cross section is calculated with a linear interpolation. In (b), (c) and (d) 
the contours for the LO cross section (obtained from HERWIG) for the SUSY background 
channel gg — > tih, gg — ► XX an d gg —>■ q,g respectively are shown. 

of constant cross section in this plot have been obtained by linear interpolation between 
the points. The NLO cross section for the SM background channel ti is 737 pb [21] • For 
the tiZ and the non-SM background channels only leading order (LO) cross section are 
used. However, the ti background is by far the largest of all the backgrounds and the 
NLO corrections to the other background channels would be small with respect to the 
tt alone. The LO cross section for the SM background channel tiZ as obtained with 
HERWIG is 439 fb. The LO cross section for the SUSY background channels depend on 
tan/? and weakly on rriA and are calculated using HERWIG with its SUSY extension |17j . 
The obtained result is shown for the SUSY background channels gg — > tih, gg — > XX an d 
gg — » q,g in figure EJo, c and d respectively. 
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4 Event Production 



Table ^ shows for each process the number of simulated events generated by HERWIG 
and processed with ATLFAST. The signal and SUSY background events were produced 
for the tan/? and via values represented by the 99 dots in figure EH- The SM backgrounds 
obviously do not depend on any MSSM parameter. 





Process 


Number Events Produced 


Signal 


gb^ 


#± - 




~4-10 5 for each (tan/?, via) 




and t 


— > fegg 


dot in figure 


SM Bkg 


99 -» 


ti 




10 8 




99 -> 


ttz 




2-10 7 


SUSY Bkg 


99 -» 


tth 




10 7 for each (tan f3, via) dot 




99 


XX 




10 7 for each (tan /3, m^) dot 




99 -» 


Q,9 




10 6 for each (tan /3, m^) dot 



Table 1: In the rightmost column the number events produced for each process, ensuring 
a negligible statistical uncertainty on the number of expected event after all cuts, are 
shown. The signal and SUSY background events were produced for the different tan/3 
and via values represented by the dots in figure EK- The SM background processes do not 
depend on MSSM parameters. 



5 Event Selection 

Four different cuts were optimized and used to enhance the signal over background ratio in 
the generated event sample. The first cut excludes events without three isolated leptons. 
This Three Lepton Cut requires: 

• Exactly three isolated leptons (£ = e or /i) with \rj\ < 2.4, with p? > 7 GeV and at 
least one of which with px > 20 GeV. 

An electron (muon) is considered as an isolated lepton if it has pr > 5 GeV (6 GeV), if 
it is within \r)\ < 2.5, if the sum of energy deposited in a hollow cone of 0.1 < AR < 0.4 
(where AR = a/ (Aip) 2 + (Ar/) 2 ) around the lepton track is less than 10 GeV and if no 
other isolated lepton nor jet is within the cone AR = 0.4 around it [THj. In order to trigger 
on the signal events the detection of high-pr leptons are required. The first cut mentioned 
above is chosen such that it meet the requirements of the ATLAS trigger system which 
has rather low pt thresholds on electrons and muons [2*2] . 
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Figure 5: Invariant mass, M«, of two isolated leptons with the same flavor and opposite 
charge is histogrammed for (a) the signal and for (b) the background channels. For the 
SUSY channels tan/3 = 10 and = 350 GeV. The histograms are normalised to the 
expected event rate for an integrated luminosity of 300 fb _1 . 



The second cut is based on the fact that for the signal two of the three isolated leptons 
come from the decay of a neutralino, x*j- The neutralino can decay to two same flavor, 
opposite charge leptons plus some undetectable particles. The invariant mass of the lepton 
pair, My, is kinematically constrained. For example when Xi is the only undetectable 
particle from x°j ( as i n figure EJd) we have [23] 

M u max = yj (m 2 .o - mfj (mf - m| ; ) /mf , (1) 

where t is the slepton involved in the decay chain x° 2£ + %\. For tan/3 = 10 

and rriA = 350 GeV we have m^o = 77.6 GeV, m^o = 131.4 GeV and m^o = 145.9 GeV 
for the MSSM point given in section H [H|. This 'gives M Umax = 50.9 GeV (68.0 GeV) 
for j = 2 (j = 3) which can be approximately seen by the first (second) edge in figure Et- 
The explanation for the tail above M« > 63 GeV is limited invariant mass resolution in 
this region. For tan/3 = 10 and uia = 350 GeV m#+ < m^o + m^+, for the MSSM 
point given in section |2l which makes a third edge (for j = 4) kinematically impossible. 
The edge corresponding to j = 3 is not very dependent on tan ft or and corresponds 
approximately to the maximum value for Mu (see figure and b). A minimum value of 
Mm = 10 GeV and a maximum value of 63 GeV is allowed. The Two Lepton Cut requires: 

• Out of the three isolated leptons at least one pair is found with two lepton of same 
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Figure 6: The sum of pt for jets that have px > 50 GeV and \r]\ < 4.5, but were not 
assigned to one of the three top jets for (a) the signal and for (b) the background channels. 
For the SUSY channels tan (3 = 10 and itla = 350 GeV. The histograms are normalised 
to the expected event rate for an integrated luminosity of 300 fb . 



flavour, opposite sign and with invariant mass in the range 10 GeV < M« < 63 GeV. 

Apart from requirements on the isolated leptons, it is also possible to impose conditions 
on the hadronic decay products of the top quark (t — » b + W) that is produced together 
with the charged Higgs boson. The Top Cut requirements are: 

• Events must have at least three jets, each with pt > 20 GeV in \r)\ < 4.5. 

• Among these, the three jets most likely to come from the top quark are selected by 
minimizing \mjjj — m t \, where rrijjj is the invariant mass of the three-jet system. It 
is required that \rrijjj — m t \ < 35 GeV. 

• Among these three top jets, the two jets most likely to come from the W boson is 
selected by minimizing \rrijj — mw\, where rrijj is the invariant mass of the two-jet 
system. It is required that \rrijj — m w \ < 15 GeV. 

The requirement that the third jet (assigning the two first to the W decay) be tagged as 
a 6-jet was also tested. The 6-tagging of jets is simulated in ATLFAST and even though 
a 6-tagging efficiency as high as 0.6 was assumed the b-Tag Cut did not increase the 
significance. 



S 




Figure 7: The significance (S/yB) after each cut is shown for tan/3 = 10 and 
wia = 350 GeV. The other MSSM parameter values are M2 = 210 GeV, /i = 135 GeV, 
mj- = 110 GeV, mf R = 210 GeV, m- g = 800 GeV, rriq = 1 TeV. An integrated luminosity 
of 300 fb" 1 is assumed. The cross sections (a) are taken from figure 0] 



The remaining jets that pass the cuts of pt > 20 GeV and \r]\ < 4.5 but which where 
not assigned to the top quark are used for the fourth cut. The q, g background contains 
many hard jets as shown in figure |Hl In view of this, the Jet Cut requires: 

• The scalar sum of the pr of all jets with px > 50 GeV, excluding jets assigned to 
the top quark, is not allowed to exceed 300 GeV. 



6 Results 

For three years of LHC operation at high luminosity (10 34 cm~ 2 s _1 ) the integrated lumi- 
nosity is about L = 300 fb" 1 . The total number of events collected during this period for 
a certain process is N to t = BR-cr-L, where the branching ratio for the signal is obtained 
from figure El and the production cross sections, a, for the signal and the backgrounds 
from figure |3J The product of the cross section and the branching ratio for gb — > tif ± 
with H — > 3/ + N can be seen in figure |H1 To achieve the expected signal rate a factor of 
2/3, for the hadronical top decay, has to be multiplied with the values in figure |H1 For the 
generation of the backgrounds the decay modes are not restricted and hence BR = 1 for 
all background processes. If the efficiency for a certain selection is e se i, then the number 
iV of events surviving this selection is iV = e se iN to t- The resulting numbers iV are given 
in table 121 for all four cuts introduced in succession for tan (3 = 10, rriA = 350 GeV, the 
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AT ■ 7 

ly signal 






Ntth 






S 

7S 


No Cut 


6068.7 


2.2-10 8 


131738 


140456 


3.310 b 


1.3-10° 


0.4 


31 Cut 


3018.3±3.3 


250938±744 


2554.0±4.1 


11.6±0.3 


241851±273 


207915±343 


3.6 


21 Cut 


2246.1±3.2 


87201±439 


129.1±0.9 


6.4±0.2 


174724±234 


82329±228 


3.8 


Top Cut 


1327.6±2.8 


12220±164 


55.2±0.6 


4.4±0.2 


4134±37 


41074±163 


5.5 


Jet Cut 


1239.7±2.7 


11995±163 


47.4±0.5 


3.2±0.2 


3092±32 


1331±30 


9.7 



Table 2: Number of events for signal and backgrounds after the successive introduction 
of the four cuts, for tan/3 = 10, itla = 350 GeV and for an integrated luminosity of 
L = 300 fb _1 . The errors quoted are derived from the Monte Carlo statistical error. In 
the last column the significance is shown. 
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values of the other MSSM parameters as specified in section and L 
significance is calculated as 

<-> _ ^signal 



300 fb" 1 . The 



(2) 



where the sum is taken over all background channels (see table [Q. In the analysis the cuts 
have been chosen so as to maximize the final significance, i.e. the significance after the 
Jet Cut for tan (3 = 10 and = 350 GeV. The progression of the signal-to-background 
rejection when imposing in sequence the four different cuts is visualized in figure [7| where 
the significance after each cut is plotted. 
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Figure 9: The contours show the selection efficiency for the signal (gb —>■ tH + ,iH~) and 
for the SUSY background channel gg — > tih, gg — > XX an d gg — > q,g in (a), (b), (c) 
and (d) respectively. The efficiencies are calculated from simulated data taken from the 
(tan /3, m^) values indicated in figure and in between the efficiency is calculated with 
linear interpolation. 

The selection efficiencies, i.e. the ratio of the number events selected by all cuts out 
of the total number events, for the signal and the SUSY backgrounds are dependent on 
the (tan/?, ttia) value. This is shown by the contours of constant selection efficiencies 
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in figure El The selection efficiencies for the SM background channels, not shown in 
figure are, for all (tan/3, m^) values, 5.4-10 -5 and 3.6-10 -4 for gg — > tt and gg — > ttZ 
respectively. 
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Figure 10: The 5-cr discovery contour for the H ± — > Xi 2X1 2 3 4 3^ + iV channel for the 
parameter set defined in section El is shown in the tan (3 vs. plane. The integrated 
luminosity is 300 fb -1 . 



The selection efficiencies are used to also calculate the significance for different 
(tan (3, vn a) values. This was done on a grid of points in the tan (3 vs. plane with 
tan f3 steps of 1 and steps of 25 GeV either directly from data at the data points 
(indicated by dots in figure HJ0 or by linear interpolation between the data points. In 
figure El the contour where the significance reaches the value 5 is shown, superimposed 
on the previously introduced discovery contours shown in figure ^ The left edge of the 
potential discovery region more or less follows the 1 ■ 10 -4 contour of branching ratio of 
~ *■ X12X1234 ~^ + shown in figure Eb- The right edge of the discovery region fol- 
low the pattern of the NLO cross section for the signal gb — > tH + , iH~ shown in figure EK 
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Figure 11: The 5-cr discovery contour for the H ± — > \\ 2X1 2 3 4 — > 3^ + N channel for the 
parameter set defined in section |21 is shown in the tan/3 vs. itla plane, superimposed over 
the corresponding result for CMS [Tl]. The integrated luminosity is 100 fb -1 . 
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7 Conclusions 



With the specific SUSY parameter set chosen, the contour of 5-cr significance for charged 
Higgs discovery through its decays to SUSY particles encloses the major part of the 
intermediate tan/5 region not covered by the charged Higgs to SM particles decays as 
shown in figure ^ The result shows a discovery region about the same size as that 
obtained for CMS ^1] using the same parameter set (the two regions are compared in 
figure ITT1 where L = 100 fb _1 ). However the regions for ATLAS and CMS differ in shape 
due to the difference in the selection criteria used, in the cross section assumptions (in the 
present analysis the NLO cross sections for gb — > tH + , iH~ were used, whereas in [Q] LO 
cross sections were used) and in the detectors used (described by ATLFAST for ATLAS 
and by CMS JET for CMS). 

As can be seen in table ED there are about ten times more ti and about three times 
more XX events than the signal events present after the last cut. Since no discriminating 
signature which helps to extract the signal from the remaining background could be 
identified a counting experiment is performed. However, in doing so, the assumption that 
the cross sections for the various background processes have been measured already is 
relied upon. 

Many more MSSM parameter sets need to be analysed before more general conclusions 
can be made, regarding the region in the parameter space where charged Higgs can be 
discovered. However, as a preliminary result, this analysis shows that for a specific SUSY 
parameter set SUSY decays of the charged Higgs bosons can be detected in ATLAS. This 
is especially interesting for the intermediate tan/? region of 4 < tan/5 < 10, where SM 
decays of the charged Higgs cannot be detected. 
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